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Abstract 

We have calculated the proton charge radius by assuming that the real proton radius is not 
unique and the radii are randomly distributed in a certain range. This is performed by averaging 
the elastic electron-proton differential cross section over the form factor cut-off. By using a dipole 
form factor and fitting the middle value of the cut-off to the low Mainz data, we found the 
lowest ^or a cut-off A = 0.8203 it 0.0003 GeV, which corresponds to a proton charge radius 

rE = 0.8333 it 0.0004 fm. The result is compatible with the recent precision measurement of the 
Lamb shift in muonic hydrogen as well as recent calculations using more sophisticated techniques. 
Our result indicates that the relative variation of the form factor cut-off should be around 21.5%. 
Based on this result we have investigated effects of the nucleon radius variation on the symmetric 
nuclear matter (SNM) and the neutron star matter (NSM) by considering the excluded volume 
effect in our calculation. The mass-radius relation of neutron star is found to be sensitive to this 
variation. The nucleon effective mass in the SNM as well as the equation of state of both the SNM 
and the NSM exhibit a similar sensitivity. 

PACS numbers: 26.60.-c, 21.65.-f,13.40.Gp, 14.20.Dh 
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I. INTRODUCTION 

The recent precise measurement of the Lamb shift in muonic hydrogen atom \^ has 
sparked a controversy, because this measurement yields a smaller proton charge radius, 
i.e. rE = 0.84184(67) fm. This radius is significantly smaller than the standard CODATA 
value {2], rE = 0.8768(69) fm, which is based on the measurements of the Lamb shift 
in electronic (conventional) hydrogen atom, as well as the results from elastic electron- 
proton scatterings. The latest precise measurement of elastic electron-proton scattering at 
MAMI, Mainz, which yields te = 0.879(5)stati4)syst.(2)modci(4)group fm [3|, clearly supports 
the CODATA value. Considerable efforts 14 -21 1 have been devoted to attack this proton 

n n 

radius problem. References |5[, for instance, propose that the off-shell form factors of 
the proton could generate large polar izability contributions to the proton structure and 
eventually could solve the problem, since the effect would only appear in the case of muonic 



Hydrogen. However, a different opinion has been put forward in Ref. [22|, in which the 
off-shell effect is shown to be not sufficiently large to reduce the discrepancy between the 
radii found in muonic and conventional Hydrogen atoms. It is interesting to note that Ref. 



22l | also concludes that the resolution of this problem must lie elsewhere, perhaps in re- 



analyses of the older experiments. Furthermore, QED is believed to be more precise than 
QCD and the techniques and methods of the Lamb shift measurement in muonic hydrogen 
as well as electron-proton scattering are beyond any doubts [23]. Therefore, it is urgent to 
reinvestigate the prevaihng methods of extracting the proton charge radius. Such an idea has 
been recently proposed by a number of research groups 2J, |25|] . Nevertheless, surprisingly 
none of them has questioned the idea of the "radius" itself. 

The radius of proton is defined in accordance with our imagination that proton has a 
spherical form. However, recent investigations have revealed that protons could deform 
from a spherical shape, like nuclei in nuclear physics. This originates from the relativistic 
motion of the spin 1/2 quarks inside the proton [26^, although there is also a claim that a 
)ure s-wave nucleon model, with / = and thus perfectly spherical, could be constructed 



271]. Meanwhile, in the liquid drop model it is also customary to assume the variation of 



the proton and neutron radii in order to explain, e.g., the polarized electric dipole moment 
in the reflection asymmetric nuclei j28|. Obviously, the definition of radius gets blurred if 
the proton were not spherical. 
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On the other hand, the fluctuating size of proton has become an important idea in 
explaining the oscillating color transparency |29|. The idea behind this fluctuating size is 
that in the proton-nucleus scattering the high energy protons that scattered at wide angles 
should be "small". However, there is also a certain process, in which the protons must 
be "large" and the amplitude of this process will increase with increasing the proton sizes. 
Using this idea the oscillating transparency found in experiment 30|, which is deflned as the 
ratio between the proton-proton scattering cross sections off the nucleus and off the proton 
at 90° as a function of energy, can be successfully reproduced [29 1. 

Based on the above experiences in this paper we propose a calculation of the proton 
charge radius by assuming that protons do not have identical radii, they vary in a certain 
range. To simplify the problem we further assume that the radii are randomly distributed 
around their average value. Practically, since the proton radius enters the cross section via 
the charge form factor, we can perform this calculation by taking the form factor cut-off as 
the corresponding variable. We believe that further corrections could enter the cross section 
formulation. However, for the present exploratory study we also believe that our assumption 
would be sufficient. 

We note that our result is in agreement with that obtained from muonic hydrogen 
A more careful measurement has been carried out at Paul Scherer Institute and the result 
has just been published 32|. It is interesting to note that the latter is still consistent with 
the previous measurement ij], indicating that the proton radius extracted from muonic 
hydrogen would hardly change. Therefore, the discrepancy between our result and the 
electronic hydrogen experiment is still outstanding and more efforts are required to alleviate 
this problem. 

Whereas a flve-percent difference in the proton radius could trigger a strong controversy 
in hadronic studies, it is quite ironic to realize that in the nuclear and neutron star matter 
investigations protons and neutrons are traditionally considered as "point particles". Ef- 
fects of the nucleon structures are only considered in the so-called excluded volume effect 
(EVE) model [33|-|39|]. In this model the total volume occupied by nucleons, i.e. Nvn, 
is subtracted from the nuclear matter volume V , so that the effective volume available for 
the nucleon motion is reduced to — Nv^, where f at is the volume of a nucleon. Of course, 
the volume of the nucleon itself decreases as the matter density increases. However, sur- 
prisingly there has been no unique deflnition of the radius in free space, i.e. at zero-density. 
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For instance, Ref. 3J| used the proton radii Vp = 0.80 fm, 0.70 fm, and 0.60 fm to study 



the EVE on the equation of state of homogeneous hadronic matter, whereas Ref. 



33| used 



0.63 fm to study the effect on the equation of state of nuclear matter. Nevertheless, 



all studies indicate that the effect is non- negligible. In fact, Ref. |3J] found that the effect 
can enlarge the range of applicability of the quark-meson-coupling model. Thus, it would 
be very interesting to study the EVE by using our knowledge obtained from the elastic 
electron-proton scattering process. 

In Sec. HI] of this paper we explain the procedure of extracting the proton charge form 
factor. Section IIIII briefly discuss the possible future experiment for refining the present 
calculation. In Sec. IIVI we investigate effects of the nucleon radius variation in the neutron 
star and symmetric nuclear matter. We will summarize and conclude our findings in Sec. El 



II. EXTRACTION OF THE PROTON RADIUS FROM ELECTRON-PROTON 
SCATTERING 



The differential cross section for elastic electron-proton scattering can be efficiently writ- 
ten in terms of the Sachs electric and magnetic form factors, Ge,p and Gm,p, as 40 1 
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where the Mott cross section. 
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J Mott EAE^ sin\e/2) 
describes the elastic scattering of point-like particle. The notation e = [1 -|- 2(1 -|- 
r) tan2(^/2)]~^ denotes the virtual photon polarization, = AEE' sin^ {6/ 2) is the square 
of the virtual photon momentum transfer, r = Q^/Arrip, irtp is the proton mass, and E [E') 
represents the electron initial (final) lab energy with scattering angle 6. 

Since we will not focus on the problems of extracting the magnetic form factor, we will use 
the phenomeno logical scaling Gm,p = fJ'pGE,p, where /ip is the proton anomalous magnetic 
moment, to simplify Eq. ([T]) to 
da 
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assuming identical charge distributions in the protons, i.e. identical electric and magnetic 
radii, where A is the corresponding cut-off. 



However, if the the proton sizes were not identical, and if we assume that they were 
randomly distributed near their middle value, then Eq. ([1]) must be averaged over all proton 
sizes, i.e., averaged over the form factor cut-off A, 
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'Ai-AA 

and 2AA represents the range of the cut-off variation around the middle value Ai. 

Experimental measurements for decades have indicated that the square root of this av- 
erage can be parameterized by means of a dipole form, 

-2 
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(6) 



where A^ = 0.71 GeV^ is often called as the standard dipole form factor, from which one 
obtains the proton electric radius by calculating the form factor slope at the real photon 
point, 

1/2 



1/2 
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dQ' 



(7) 



At this stage it is important to note that after the operation of modern continuous beam 
accelerators, such as MAMI in Mainz and CEBAF at the Jefferson Laboratory, significant 
deviation from the standard dipole form factor has been observed. To account for this 



deviation a number o 
of the dipole form 



43 



new fits and models has been proposed. This includes modifications 



factor 



42 



4J], as well as introduction of more physical ingredients in the form 



45l | . There seems to be no need to keep the original dipole form to fit both electric 



and magnetic form factors, especially for a global fit to all data, since the standard dipole is 
considered as just a phenomeno logical approximation. Furthermore, the choice of the dipole 
form also seems to be trivial. 

However, in the non-relativistic limit as well as in the Breit frame the dipole form factor 
is related to an exponentially decaying charge distribution via the Fourier transform. The 
exponentially decaying behavior is found in most natural phenomena, from radioactive decay 
rate to atmospheric pressure on earth. Some phenomena in social science also exhibit this 
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behavior. Thus, we beheve that in our present case a dipole form factor looks more natural 
and a deviation from such a natural phenomenon requires a rigorous physical concept. 

Furthermore, it is also important to emphasize here that the determination of the slope 
at = given by Eq. (j?!) requires very good knowledge of Ge,p{Q'^) at very low Q^. As 
a consequence, the closer we can approach the real photon point experimentally, the more 
reliable we can determine the proton charge radius. Therefore, the latest and accurate 
measurement of the elastic electron-proton scattering with low energy and low at MAMI 
[st provides very suitable experimental data for our present discussion. 

We begin with Eq. (E]), which implies that the experimentally observed form factor is in 
fact an average to the genuine form factor, that one actually should use in Eq. ([7]) in order 
to get the real proton radius. As a consequence, the extracted proton radius in this way 
should be considered as an averaged radius. 

To investigate the effect of averaging the form factor given by Eq. ([5]), let us use the 
standard dipole form factor to calculate Ai)) in Eq. ([5]). Note that if we use a 

dipole form, the magnitude of the relative variations of both radius r and cut-off A are equal, 
i.e. 

A A A „ 

(8) 

provided that the variations are not extremely large. 

The result is shown in Fig. [H where we compare our calculations obtained with the cut- 
off variations AA from zero up to 40% of its standard value with experimental data. Note 
that for the sake of simplicity we use the latest result from Mainz experiment jsj], which 
provides the latest and most accurate data in the low region, and the result of extraction 
rom the world electron-proton scattering data with two-photon exchange effects included 



AA 




Ar 


A 




r 



42( 1 . which represents the previous measurements. We use the result obtained from the 



standard Rosenbluth separation technique for the Mainz data, in order to reduce the model 



dependency of the data. The Ge,p data extracted in Ref. [42j are of course model dependent. 
However, in this paper they are only used for the purpose of comparison and not included 
in the fitting database as described in the following discussion. 

It is obvious from Fig. [1] that the result shows a variance to the standard dipole form. 
For > 0.25 GeV^ we observe that the averaged form factors are larger than the standard 
dipole one (i.e. AA/A = 0); increasing the relative variation will increase the form factor. 



6 




(GeV2) 



FIG. 1: (Color online) Square root of the averaged calculated from Eq. ([5]) for different 

values of the relative variation AA/Ai (shown in the figure) compared with experimental data. 
Experimental data are taken from Refs. 3|] (open circles) and [42] (solid squares). 



However, for < 0.25 GeV^ we observe a different behavior, i.e. the form factor decreases 
as the variation increases. Since the decrease is relatively small, it is almost invisible in 
Fig. [U Therefore, in Fig. |2] we increase the resolution of the ((^^^(Q^, Ai))^/^ axis by 
limiting < 0.25 GeV^, where we can clearly see the effect of variation, i.e. for the relative 
variation of 40% the agreement with experimental data is almost perfect. 

Although the agreement of the solid curve with experimental data in Fig. [2] could be 
fortuitous, the most important message is that the standard dipole form factor is still valid 
at low Q^, provided that the corresponding cut-off must be averaged with relative variation 
AA/Ai = 40%. We believe that this is crucial because the extraction of proton charge 
radius is always plagued with many complicated corrections, especially at high regime, 
as discussed above. In view of this, in what follows, we will only use the MAMI data and 
limit the only up to 0.25 GeV^. 

It is apparent from Eq. ([5]) and Fig. [2] that for each value of AA we can optimize the 
form factor cut-off Ai in order to further improve the agreement of our calculation with 
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(GeV2) 

FIG. 2: (Color online) As in Fig. [H but limited for < 0.25 GeV^. The results obtained for two 
different values of AA/Ai are compared with the standard dipole form factor and experimental 
data. Note the linear scale for axis. 

experimental data. For this purpose we can calculate the standard 'X^ /N which measures 
the agreement of our calculation with experimental data. The result as functions of the 
proton radius (translated from Ai) and AA/Ai is displayed in Fig. [3l It is obvious that the 

/N has only one minimum located by the intersection of the two dashed lines. 

To locate this minimum accurately we fit the value of Ai by using the CERN-MINUIT 
code and scan the relative variation AA/Ai from 0% to 50% with 2.5% step, simultaneously, 
where the value of Ai is allowed to vary between 0.80 and 0.90. Although this choice seems 
to be arbitrary, in our fits we found that the Ai value never reaches both upper and lower 
limits. For AA/Ai = 0% (50%) the cut-off value is obtained to be 0.8142 (0.8493) GeV, 
which corresponds to the proton charge radius of 0.8396 (0.8048) fm. 

Having finished the scanning process we observe that the obtained 'X^ /N forms a parabola 
with the minimum value at AA/Ai = 21.5% and Ai = 0.8203 GeV. This corresponds to 
the relative variation of proton radius Ate/te ~ 21.5%, from Eq. ([8]), and the proton 
radius of r = 0.8333 fm, from Eq. ([71). The complete result of this scanning process is 
depicted in the left panel of Fig. HI To increase the accuracy, we have refined the AA/Ai 
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FIG. 3: (Color online) (a) The obtained /N as functions of the proton charge radius r and 
A = AA/Ai. (b) The projection of /N on the r — A plane for /N < 16. Intersection of the 
vertical and horizontal dashed lines locates the minimum position of x? /N , which is accurately 
shown in Fig. 



step to 0.005, which is equivalent to Ar^; = ±0.0003 fm, in the vicinity of the minimum. 
Therefore, our calculation would produce the best agreement with experimental data if we 
used te = 0.8333 ± 0.0004 fm, where we have added the error bar coming from the fitting 
process (Ar obtained from the MINUIT package). The present result is very interesting 
)ecause it corroborates most of the latest findings that exploits more sophisticated techniques 



24 



25|. 



We have also performed the above procedure to find the proton magnetic radius. It is 
well known that the experimental data in this case are notoriously inaccurate, especially at 

^0. As a consequence, we did not use the four lowest data points from Ref. 31| . 
because most of them cannot be renormalized to yUp at the real photon point. For comparison 
with the charge radius, we display the result in the right panel of Fig. |H Obviously, the 
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FIG. 4: (Color online) The /N as a function of the obtained proton (a) charge and (b) magnetic 
radii. The lowest and highest values of relative variation in the form factor cut-off, A = AA/Ai, 
are indicated in the figure. The vertical lines indicate the minimum positions of /N . The 
corresponding A values are also shown. 



trend is different as in the case of the charge radius. In the case of the charge form factor, 
increasing the relative variation is required to decrease the magnitude of Ge,p in order to 
reproduce experimental data (see Fig. [2]). In contrast to this, the relative variation of the 
radius is required to increase the magnitude of Gm,p, in order to reproduce the data. 

As shown in Fig. IHwe obtain vm = 0.8103 ± 0.0004 fm, which is smaller than the result 
extracted from the dispersion relation, i.e. 0.84^qq2 ['^'^l- Nevertheless, our magnetic 
radius is much larger than that obtained from direct extraction of the Mainz data, i.e. 
0.777 ± 0.013 fm [3,]. However, if the Friedrich-Walcher parameterization 43|] was used in 
the latter, the magnetic radius would increase to 0.807 ± 0.02 fm 3lj, which is apparently 
in good agreement with our finding. We believe that the less accurate magnetic form factor 
extracted from the electron-proton scattering could be the origin of the large variance in the 
extracted magnetic radii found in the literature. 
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III. EXPECTED FUTURE ELECTRON-PROTON SCATTERING EXPERI- 
MENTS 



Since the mathematical formula of the proton charge form factor is in principle not known, 
determination of the proton radius using Eq. ([7]) requires very good knowledge of the proton 
form factor to a very low region. Thus, the real challenge for future experiments is 
to extend the current experimental data to this kinematics. The situation is exhibited in 
Fig. [5l where we compare the result for the 21.5% proton radius variation obtained in the 
previous section and various available parameterizations with experimental data js, 42|. It 
is obvious from this figure that the Mainz data tend to deviate from our present result, 
whereas, surprisingly, the data extracted in Ref. 42| show a very good agreement with our 



calculation. We note that Ref. 



42| used polynomial expansion to parameterize the form 



factor during the extraction. Therefore, we believe that the agreement with the present 
calculation as exhibited in Fig. [5] could not be a coincidence. 

As in the large case, in the very low region it is also obvious that the standard dipole 
form factor is substantially larger than our calculation. The result of our calculation is very 
close to the result obtained from the Friedrich-Walcher model 43|. Since all four models 
shown in Fig. [5] yield significantly different proton charge radii, it is of course important 
to refine the experimental measurement at this kinematics. Theoretically, this is possible 
because Eq. ([1]) explicitly shows that at low contribution of Ge-p is dominant. Thus, at 
low measurement of Ge,p should be more accurate than that of Gm,p- 

However, from the experimental side this could be a daunting task. It should be remem- 
bered that measurements of electron-proton scattering cross section for ^ 0.05 GeV^ 
were already at forward angles 3l|- Below this point, presumably one has to use other 
methods. One possible choice proposed at MAMI is the use of the initial state radiation, 
i.e. radiation emitted by electron before it is scattered by the proton, which could provide 



measurement of form factors down to = 0.0001 GeV^ 
could work, the proton radius would be severely constrained. 



3l[ |. Obviously, if this method 
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FIG. 5: (Color online) Square root of the averaged ^ obtained from different calculations 
compared with experimental data for very low Q^. Notation for experimental data is as in Fig. [TJ 
Result obtained in the present work is given by the solid curve. Parameters for the double-dipole 



and Friedrich-Walcher form factors are taken from Ref. 3l| . 



IV. EFFECTS OF THE NUCLEON RADIUS VARIATION ON THE SYMMET- 
RIC NUCLEAR MATTER AND THE NEUTRON STAR MATTER 

In the relativistic mean field (RMF) model the Lagrangian density of nucleons consists 
of four terms, i.e. the free nucleon, free meson, interaction between nucleons via meson 
exchange, and meson self interaction terms. If we assume that electrons and muons are 
point particles, whereas nucleons have structures with a radius r^, then according to the 
RMF model the energy density of matter consisting of the nucleons and the leptons is given 



by [39 



4ll 



e = A(eJ + e^) + + + eM(w, a, p) 

+ 9uj(^oipp + Pn) + -gpboiPp - Pn), (9) 

where g^, and gp are the couplings for a and p mesons, respectively, e^/ is the total 
energy density of the meson, while cr, ujq and 6o are the cr, u and p fields, respectively. 
Furthermore, in Eq. ([9]) we have 

e-=-^[d'k{e + m*'f'eik-kF), z = p,n,e,/i, (10) 
(27r) J 
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where for leptons the effective mass is m* = rrii and for nucleons m*= rrii — g^a. 
The nucleon and scalar densities read 



Pi = Api , (11) 
Ps,i = Aps,i , (12) 

where pi and ps^i are the ith nucleon and scalar densities, assuming the nucleon is a point 
particle. The normalization constant A is given by 

A = , (13) 

l + VpPp + VnPn ^ ' 

with Vp and Vn the proton and neutron volumes, respectively. To simplify the present 
calculation we assume that 

Vp = Vn = Vi, = ^vrr^, (14) 

where Vjq and r^v are the volume and radius of the nucleon, respectively. 
From Eq. (|9]) we can derive the matter pressure, 

with e = t/ p. Furthermore, the chemical potential for the ith nucleon can be obtained from 

Pi = Ep ^ + ViPl + g^ujQ + ai]^gphQ, (16) 
with ctj equals +1 (—1) for proton (neutron), E*p^ = {kp^ + m*^)"^^^ and 

3 *4 , / ^F,i + E*p i 



Different from the quark meson coupling (QMC) model 3J, |36|, where the dependence 
of ttv on the matter density can be directly obtained from the model, in the RMF approach 
the dependence cannot be easily predicted. Therefore, in the present study we choose a 
phenomenological form for the nucleon radius, which is given by 

rN{p)=rN{^){l + Hj] \ ' (18) 
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FIG. 6: (Color online) Comparison between proton bag radii as a function of the ratio between 



nucleon and nuclear saturation densities obtained from the QMC models 
two different /3 values. 



Mi and Eq. ([HD with 



where p = pp + pn, po is the value of p at saturation point, and r^iO) is the proton radius 
in vacuum (zero density), determined from Eq. ([7]). At first glance, the choice seems to 
be trivial. However, it is actually selected to fulfill the causality constraint. Furthermore, 
the formula given in Eq. (fT8|) is more convenient for the present purpose, rather than the 
exponential one, because in the framework of the presently used RMF model we found that 
the required radius must slowly fall off as a function of density. Otherwise, the predicted 
neutron star mass would violently overshoot the mass of the SRJ164-2230 pulsar, which is 



believed to be the heaviest observed neutron star |46j. We also observed that, for selected 
value of /3, Eq. (fT8|) can be adjusted to mimic the result of the QMC models |34] in a certain 
range of density, thus providing a good check of our result. 

The ratio between r^lp) and r^iO) is exhibited in Fig. |6l where we compare the results 
obtained from two f3 values, i.e. (3 = 0.0005 and 0.01, with those obtained from the QMC 
calculation using different values of bag constants B (see Ref. SJ] for explanation). It is 
apparent from this figure that the difference between the calculated radii increases as the 
density increases. Nevertheless, the distributions of the radii obtained from the QMC model 
are still bounded within the difference of the two chosen /3 values in Eq. ( !T8|) . i.e. the solid 
and dash-dotted lines in Fig. |6l 
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FIG. 7: (Color online) The neutron star mass as a function of its radius for four different nucleon 
radius assumptions. Panels (a) and (b) show the calculation with /3 = 0.0005 and (3 = 0.01 in 
Eq. (I18p . respectively. The gray horizontal bands show the mass of the SRJ164-2230 pulsar, which 



is believed to be the heaviest observed neutron star 
given in Fig. [HI 



46l |. Notation of the curves in both panels is 



For the RMF model we use the parameter set obtained by the lUFSU collaboration |4j 
In calculating the equation of state (EOS) of the neutron star matter (NSM) we use the 
neutrality and /3-stability conditions. They are required in calculating the Fermi momentum 
of each particle in the neutron star core. 

The neutron star mass as a function of its radius can be obtained by solving the Tolman- 
Oppenheimer-Volkoff (TOV) equation with different particle densities in the neutron star 



core. 
et al. 



n order to describe the "outer crust" region we have used the EOS given by Riister 



48| . The EOS for the "inner crust" region is obtained from an extrapolation of the 



EOS of the core and "outer crust" by making use of the polytrophic energy-pressure density 
approximation. The "crust" region is described by using the RMF model with and without 
the EVE. 

Since in the nuclear and neutron star matter a direct comparison between model calcu- 
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lations and precise experimental data, as in the case of electron-proton scattering in the 
previous section, is beyond our imagination at present, in the followings we will not calcu- 
late the effect of averaging the nucleon radius on the possible observables. Furthermore, as 
already obtained in the previous section, the effect of this averaging process is a relatively 
tiny shift from the original value. Therefore, we believe that at this stage it is sufficient 
to investigate the effect on the conventional observables by using a ±20% variation of the 
nucleon original radius. Note that this variation enters our calculation through Eqs. f|T3|) 
and ([HD. 

Figure [7| shows the sensitivity of the neuron star mass-radius relation to the variation of 
the nucleon radius as well as to the dependence of the nucleon radius on the matter density 
[P in Eq. flTSl) ]. in the framework of RMF models. The radius of the neutron star depends 
on the value of nucleon radius, whereas the maximum mass of the neutron star is controlled 
by the dependence of the nucleon radius on the matter density. Therefore, the heaviest 
observed neutron star mass PSRJ1614-2230 46|] shown by the horizontal gray lines in Fig. [7] 
yields a significant suppression of the nucleon radius at very high density. This result is 
interesting, because it opens the possibility of hyperon existence in a neutron star by using 
the hyperon vector couplings obtained from SU(6) symmetry 49 1. 

With regard to the radius of the canonical neutron star (I.4M0), which is constrained 



between 10.4 km and 12.9 km 



50| , our present result should be carefully interpreted, because 



in this calculation we have used the lUFSU parameter set, which was fitted to the finite nuclei 
data by assuming point particle approximation for the nucleon. Therefore, the present result 
cannot be quantitatively compared with those obtained with other constraints. Nevertheless, 
Fig. [7] indicates that the present result could become a stringent constraint to the nucleon 
radius at high density, once a consistent EOS with EVE were available. 

Since the result obtained by using /3 = 0.0005 substantially overshoots the PSRJ1614- 
2230 constraint, as shown in Fig. [3, in the following discussion we will only use (3 = 0.01. 
The result for the effective nucleon mass in the case of symmetric nuclear matter (SNM) 
is shown in Fig. [HI where we compare the calculated masses obtained by assuming point 
particle approximation [r7v(0) = 0] and finite nucleon radii [r7v(0) 7^ 0]. 

From Fig. [8] it is apparent that at high densities the nucleon effective mass increases with 
increasing the nucleon radius. In view of the instability against the particle-hole excitation at 
high densities due to the density fluctuation [5l|, a sufficiently large effective mass predicted 
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FIG. 8: (Color online) Effective nucleon mass as a function of the ratio between nucleon and nuclear 
saturation densities. The results are obtained with different values of the nucleon radius rAr(O) as 
indicated in the figure (in the unit of fm). 




FIG. 9: (Color online) Equation of states of the symmetric nuclear matter obtained from calcu- 
lations with different values of the nucleon radius as a function of (a) energy and (b) density. 
Notation of the curves is as in Fig. [51 



by the RMF model has an obvious advantage. 

Figures M and [10] exhibit the non-zero nucleon radius effect on the EOS of SNM and 
NSM, respectively. The solid circles in the NSM EOS of Fig. [10] indicate the positions of the 
NS center pressures and center energy densities of maximum mass. At moderate densities, 
which correspond to e < 300, we observe that the EOS becomes stiffer as the nucleon radius 
increases. Beyond this range, the EOS tends to be softer. From Fig. [TTl it is obvious that 
the pressure of the NS center with maximum mass of r^lO) = 1 fm (blue dashed line) is 
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FIG. 10: (Color online) As in Fig. [9^, but for the neutron star matter at (a) lower and (b) higher 
energies. The dash-dot-dotted line in panel (a) shows the causality constraint, the solid circles 
indicate the center pressures and energy densities in the neutron star with maximum mass. 

approximately 32 MeV fm~^, which corresponds to the energy density of about 320 MeV 
fm~^. At this point the EOS obtained with rAr(O) = 1 fm is stiffest compared to other 
cases [panel (a) of Fig. 10]. Since in obtaining the NS mass we should integrate the TOV 
equation using the corresponding EOS as input from the NS center pressure up to zero, 
information based solely on the soft EOS at high densities of riv(O) = 1 fm is insufficient 
for a complete understanding of the NS maximum mass. For other non zero nucleon radius 
cases, the situation is similar. 

Obviously, decreasing the nucleon radius will decrease the pressure in the region of e < 
300 MeV fm~^. However, decreasing the nucleon radius will simultaneously increase the NS 
center pressure. Thus, additional contribution from the center pressure up to 300 MeV fm~^ 
will only slightly increase the NS mass. This is due to the fact that the corresponding EOS 
at high densities are relatively soft. On the other hand, in the case of zero nucleon radius 
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FIG. 11: (Color online) Neutron star mass as a function of its center pressure for different nucleon 
radii obtained with /? = 0.01. Solid circles indicate the maximum masses with the corresponding 
center pressures. 



but using the same RMF parameter set, the contribution from 300 MeV fm~^ up to zero is 
very small but the contribution from high densities is dominant because the corresponding 
EOS is stiffer than that of the nonzero nucleon radius. The very small contribution in the 
region < 300 MeV fm~^ and the relatively large center pressure are typical for point-particle 
RMF models. In this case, the strong correlation between the NS maximum mass and the 
EOS stiffness at high densities is very obvious. Therefore, in the point particle case we 
only need to consider the EOS at high densities in investigating the NS maximum mass 
behavior. In addition, the high densities EOS become significantly stiffer as /3 decreases. 
As a consequence, decreasing the /3 value will increase the maximum mass of NS. Thus, the 
increasing of the NS mass depends sensitively on the nucleon radius in free space. These 
phenomena explain the increase of the predicted mass and radius of the neutron star with 
increasing the nucleon radius, as shown in Fig. [71 

In Fig. [12] we display the finite nucleon radius effect on the fraction of the matter con- 
stituents in the neutron star. It is obvious from this figure that increasing the nucleon radius 
will increase the number of existing charge particles at high density. This result has a seri- 
ous consequence on a number of neutron star properties, such as the neutron star stability, 
neutrino transport in the neutron star, as well as the cooling process of a neutron star. 
Unfortunately, a more detailed and quantitative analysis of the EVE on neutron star should 
wait for a more consistent EOS, which includes the EVE in the calculation. 
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FIG. 12: (Color online) Proton, neutron, electron, and muon fractions in the neutron star matter 
as a function of the ratio between nucleon and nuclear saturation densities. Notation of the curves 
is as in Fig. [HI 



V. SUMMARY AND CONCLUSION 

We have investigated effects of the proton radius variation on the extraction of the proton 
charge form factor. To achieve the best agreement with experimental data we have averaged 
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a dipole form factor over the corresponding cut-off with an upper (lower) integration hmit of 
+21.5% (—21.5%) from its middle value. The extracted proton charge radius is found to be 
smaller than that obtained using the traditional standard dipole fit, but is in good agreement 
with those obtained from a recent measurement of the Lamb shift in muonic hydrogen atom 
as well as from the dispersion relation. The extracted proton magnetic radius is smaller than 
the result of the dispersion relation, but in agreement with the direct extraction by making 
use of the Friedrich-Walcher form factor. Nevertheless, since the magnetic form factor is less 
accurate, the extraction of magnetic radius is also less reliable as compared to the result of 
the charge radius. 

We have also investigated effects of the nucleon radius variation on the SNM and the 
NSM. To this end, the nucleon radius dependence on the matter density is described by a 
simple phenomenological form and four assumptions of the nucleon radius at zero density 
are considered in the calculation, i.e. fm (point particle approximation), 0.833 fm (original 
radius), as well as 0.667 fm and 1.000 fm (±20% modifications of the original radius). We 
found that the relation between the mass and radius of a neutron star is very sensitive to the 
radius of the nucleon. A similar result is also observed in the case of the effective nucleon 
mass of the SNM as well as the EOS of both NSM and SNM. However, a more quantitative 
conclusion could be drawn only after a more consistent EOS, with the EVE considered, had 
been available. 
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